Effects of Axial Vorticity in Elongated Mixtures of Bose-Einstein Condensates 
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We consider a meniscus between rotating and nonrotating species in the Bose-Einstein condensate 
(BEC) with repulsive inter-atomic interactions, confined to a pipe-shaped trap. In this setting, we 
derive a system of coupled one-dimensional (ID) nonpolynomial Schrodinger equations (NPSEs) for 
two mean-field wave functions. Using these equations, we analyze the phase separation/mixing in 
the pipe with periodic axial boundary conditions, i.e. in a toroidal configuration. We find that the 
onset of the mixing, in the form of suction, i.e., filling the empty core in the vortical component by 
its nonrotating counterpart, crucially depends on the vorticity of the first component, and on the 
strengths of the inter-atomic interactions. 

PACS numbers: 03.75.-b, 03.75.Lm, 64.75.-g 



Since the creation of vortices in Bose-Einstein conden- 
sates (BECs) [11,121, this topic has been a subject of many 
experimental and theoretical works, as reviewed in Refs. 

In particular, much attention has been drawn to vor- 
tices in mixtures of two BEC species; in fact, the first vor- 
tices were created in a two-component setting [![, and a 
theoretical analysis of that setting was developed too [1] . 
In this connection, a situation of straightforward interest 
is the interaction of rotating and nonrotating immiscible 
BEC species. It is natural to expect that the nonrotat- 
ing component may fill the hollow vortical core(s) in the 
rotating one. As shown experimentally, vortices and vor- 
tex lattices with empty and filled cores feature a great 
difference in their structure 0, Q ■ Very recently it has 
been also predicted the possibility to create vortex with 
arbitrary topological charge by phase engineering [6(. 

Matter-wave vortices can be created not only in large- 
aspect ratio settings, but also in narrow cigar-shaped 
traps ("pipes"), which help to stabilize them In the 
pipe geometry, it is interesting too to consider the inter- 
action of rotating and nonrotating BEC species, which is 
the subject of the present work. In particular, a natural 
issue in this case is the effect of "suction" , i.e., onset of ef- 
fective mixing between two nominally immiscible species 
by pushing the nonrotating component into the empty 
core in the vorticity-carrying one. In other words, the 
suction implies indefinite stretching of the meniscus sep- 
arating the immiscible species which originally fill two 
halves of the tube. 

In the mean-field approximation, the starting point of 
the analysis is a system of coupled Gross-Pitaevskii equa- 
tions (GPEs) for macroscopic wave functions, tpi and ip2, 
of the two BEC species confined in the tight cylindrical 
trap. In the scaled form, the equations are 
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where z is the axial coordinate, while x and y are the 
transverse ones. In these equations, lengths are mea- 
sured in units of a± — \/h/ (miili), the energy in units 
of hfli, and time in units of l/fii, where and Q.1,2 
are masses and transverse trapping frequencies of the two 
species, while the relative parameters in Eqs. ([1]) and 
([2]) are m = 1112/ mi and Q = f^/f^i- The interaction 
strengths in the equations are expressed in terms of the 
s-wave inter-atomic scattering lengths, 31,2 = 2ai i 2 / a± , 

gi2 = 2ai2/<2x , where aj is the scattering length in the 
j-th species, while a\2 is the scattering length for atoms 
belonging to the different species. In this work, we con- 
sider the most relevant situation with the repulsion be- 
tween atoms belonging to the same and different species, 
i.e., we take 51,32,512 > 0. 

We will solve equations (Q]) and @, as well as effec- 
tive ID equations to be derived from them, with periodic 
boundary conditions in the axial direction (z), which as- 
sumed that the pipe is closed into a torus of a large ra- 
dius. Toroidal traps of the magnetic type are currently 
available to the experiment [8[ , and the BEC dynamics in 
the toroidal geometry was studied theoretically in several 
works [?|. Actually, a variety of the toroidal configura- 
tion is a skyrmion pattern predicted in a two-component 
BEC, in which one component, which carries the vortic- 
ity, is effectiviely trapped in a doughnut region created 
by the other component, the entire configuration being 
stable [nj. 

Equations |T]) and ^ conserve the norms of the two 
wave functions, i.e., numbers of atoms in the species, 
(1) f l^i^Cr, t)\ 2 dr = Ni,2- The conserved energy (Hamil- 
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It is well known that, in the case of repulsion between 
atoms, the two species in the free space are miscible un- 
der the condition gf 2 / (3132) < 1, and miscible in the op- 
posite case fTlj | . The pressure induced by the transverse 
confinement changes the situation, pushing the system 
towards stronger miscibility (see, e.g., Ref. fl2j]). 

We aim to study a different effect, viz., a shift of the 
miscibility threshold induced by vorticity imparted to 
one of the species. First, we derive a system of effective 
one-dimensional (ID) nonpolynomial Schrodinger equa- 
tions (NPSEs) for the rotating and nonrotating species 
trapped in the pipe. The derivation extends the earlier 
developed analysis of the single-component [l3| and two- 
component [l4| settings, as well as the single-component 
case with the vorticity [TQ| . Then, the coupled NPSEs are 
used to predict, in an analytical form, the threshold for 
the transition to the immiscibility in the mixture includ- 
ing the rotating and nonrotating components. Finally, 
the prediction is compared to direct numerical solutions 
of the 3D equations, (HJ and which demonstrates a 
high accuracy of the analytical miscibility condition. Ac- 
tual mixed and phase-separated configurations, and the 
onset of suction, are presented using numerical solutions 
of the 3D equations. 

Coupled GPEs (fTJ) and can be derived from the 
action functional, 
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with E taken as per Eq. (|3|). To derive effective ID 
equations in the axial direction (z), we generalize the 3D 
ansatz proposed for nonrotating configurations in Ref. 
[Til ], assuming that the vorticity, which is quantified by 
integer positive "spin" S, is imparted to the first species, 
while the second species has no vorticity: 
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where ai^(z, t) are axially nonuniform transverse widths 
of the two components, and fi${z,t) are the respective 
axial wave functions. We substitute this ansatz in Eq. 
P| and perform the integration over x and y, neglecting 
the derivatives of a\ t 2(z,t). In this way, we arrive at an 



effective action functional. Extremizing this functional 
with respect to f*(z,t) and f 2 {z,t) (the asterisk stands 
for the complex conjugation), we derive the coupled time- 
dependent NPSEs, 
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In addition, two relations are generated by varying the 
effective action functional with respect to a\^{z,t), 
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In the absence of vorticity in the first species, 5 = 0, Eqs. 
0) , © , ([9]) , (HO]) reduce to those recently derived in Ref. 
[lj] for the nonrotating two-component mixture. 

In the case of gi2 = the two components decouple 
and we obtain the NPSE for the single-component BEC 
under transverse harmonic confinement, carrying axial 
vorticity S [l5|. With S = 0, the equation reduces to 
the NPSE originally derived in Ref. The single- 

component model with S =/= was recently studied in 
Ref. [H. 

To analyze how the well-known phase-separation con- 
dition in the free space is modified in the pipe-shaped 
toroidal trap (i.e., as said above, we assume periodic 
boundary conditions in the direction of z), we consider 
the following initial expression for a weak perturbation 
of the axially uniform configuration (i.e., a fully mixed 
one), /i,2(z) = y/ n i,2 [1 ± a cos (z/R)], where n\ and 112 
are densities of the two components in the unperturbed 
state, a and (3 the perturbation amplitudes, and R the 
radius of the torus. Inserting these functions into expres- 
sion we find a correction to the effective energy at 
the second order in a and /?, 
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where widths <j\ and a 2 are taken as per Eqs. ^ and 
([TO]) with constant densities ri\ and n 2 . 

The uniform configuration, with a = /? = 0, is stable 
against the phase-separating axial density modulations, 
i.e., the binary BEC is miscible, if the respective energy 

d 2 E 2 d 2 E 2 



curvature, 



W 1 



\ dad/3 J 



, is positive at a = /3 = 0. 
This condition was recently applied to the investigation 
of the stability of various mixtures, viz., 3D Fermi- Fermi 
[l6| and Bose-Fermi- Fermi ones in terms of the BCS-BEC 



crossover 



17| . and ID Bose-Fermi- Fermi mixture in the 



Tonks-Girardeau regime [l8| . 

The critical curve of the phase separation is determined 
by equating the energy curvature to zero, which yields 
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From this condition it follows that, even for g\ = g 2 = 
(no intra-species repulsion, while the inter-species repul- 
sion is in the action) , the mixed state exists if the value of 
the coefficient gi 2 accounting for the repulsion between 
the species, is below a critical value, i.e. if: 
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This is clearly a finite-size effect, since expression (fT2|) 
vanishes at R — > oo . 

If the density is small enough, relations ([9]) and (fTO]) 
may be approximated by <j\ » 1 and a 2 ~ l/(m 1 / 4 r2 1 / 2 ), 
and Eq. (|llj) yields a simple explicit result, in the limit 
of R — > oo, for the onset of the "suction" effect (uniform 
density of the two species along z): 
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In the absence of axial vorticity (S = 0), and for the 
two species with the same mass (m = 1) and common 
strength of the transverse confinement (17 = 1), Eq. (fT3]) 
coincides with its counterpart for the free (3D) space, 

Svxl (5152) = 1. 

We have tested condition (fT3"|) against direct numerical 
solutions of the full 3D equations, ([I]) and (|2|). Results 
of the comparison are displayed in Fig. [TJ which shows 
that the analytical criterion (fTB")) for phase separation is 
very accurate, at least up to g\N = 10. 

We report also typical demixed (phase-separated) and 
mixed states, obtained as direct numerical solutions of 
Eqs. 0, ©, ©, HUD with periodic boundary condi- 
tions along z, to model the torus of a large radius. If 




FIG. 1: Critical phase-separation curves for the bosonic mix- 
ture, with m = fi = 1, N\ = N2 = N, and g\ = gi- Lines: the 
prediction of Eq. (|13[) : symbols: results of numerical solutions 
of Eqs. (H| and ©. 
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FIG. 2: (Color online). The axial, Pj(z), and radial, Pj(r), 
densities [see Eqs. (|14[1 and (|15p ] of the two components of 
the binary BEC in the quasi- ID toroidal trap. The continu- 
ous and dashed curves pertain, severally, to the first (carrying 
vorticity S) and second (zero-vorticity) components. The ax- 
ial coordinate takes values —20 < z < +20, with periodic 
boundary conditions. Here, m = Q = 1, Ni = N2 = N, 
giN = g 2 N = 1, g 12 N = 5. 



the radius is much larger than the transverse width a±, 
effects of the curvature in the axial direction may be ne- 
glected [l||. We solved the equations numerically by 
means of a finite-difference Crank-Nicholson predictor- 
corrector method with the cylindrical symmetry (details 
of the method were given in Ref. [13] )■ To generate the 
ground-state wave functions ipi,2{r, z), with j = 1, 2, the 
imaginary-time integration was used. 

In Fig. [2] we plot typical examples of the axial and 
radial densities, 
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FIG. 3: (Color online). The axial densities Pj(z) of the two 
components of the binary BEC in the quasi-lD toroidal trap 
for two values of the axial length, L = 40 and L = 400. The 
continuous and dashed curves represent the first and second 
components. The parameters are m = Q = 1, Ni = N2 = N, 
giN = g 2 N = 1, g 12 N = 30, S = 0. 

/+00 
dz\ih, 2 (r,z)\ 2 , j = l,2. (15) 
-00 

Figure [5] displays the demixing (phase separation) for 
5 = and 5 = 1, while no separation is observed in 
the axial direction for 5 = 2. In fact, this means that 
the empty core induced by the vorticity with 5 = 2 is 
wide enough to produce the "suction" effect, while 5 = 1 
is not sufficient for that. 

To further illustrate the situation, in Fig. [3] we plot 
the axial density profiles of the two components of the 
binary BEC in the regime of the phase separation for 



two different values of the axial length, L = 2ttR = 40, 
and L = 400, i.e., "short" and "long" tori, respectively. 
The figure shows that the density profiles get flatter and 
interfaces steeper with the increase of L. In fact, for 
very large L the contribution (generated by the gradient 
terms) to the total energy §5§ from to the interfaces be- 
come negligible in comparison with the bulk terms (the 
Thomas- Fermi approximation), the corresponding axial 
profiles looking as step functions. In this limit, the onset 
of phase separation is accurately described by Eq. (fT3)) 

In conclusion, we have considered dynamical states in 
the binary BEC formed by two species with repulsion be- 
tween atoms, in the case when one species is prepared in 
a vortical state, with vorticity 5 = 1 or 2, while the other 
has zero vorticity. It is assumed that the condensate was 
loaded in a quasi-lD toroidal trap, i.e., the correspond- 
ing equations were solved with periodic boundary condi- 
tions in the axial direction. In this setting, we have de- 
rived coupled ID nonpolynomial Schrodinger equations 
(NPSEs) for the mean-field wave functions of the two 
components. Using these equations, we have found the 
phase-separation threshold, as a condition for the onset 
of the modulational instability of axially uniform states. 
The predictions produced by the NPSEs were compared 
to results of direct simulations of the underlying 3D equa- 
tions, demonstrating a very good agreement. Further, 
stable 3D states with the species mixed and separated in 
the axial direction were found, the transition to the effec- 
tive mixing being accounted for by the suction, i.e., filling 
the empty core in the vortical component by the nonro- 
tating one. The onset of the suction depends on the size 
of the vorticity in the first component, and on the rela- 
tive strength of the inter-species repulsion in comparison 
with the intrinsic repulsion in each component. 
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